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O ' Abstract 

i> : 

. With the help of the factorizing F-matrix, the scalar products of the Uq{gl{l\l)) 
free fermion model are represented by determinants. By means of these results, we 

l/-^ ' obtain the determinant representations of correlation functions of the model. 
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I Introduction 



The computation of correlation functions is one of the challenging problem in the theory 
of quantum integrable lattice models [1, 2]. In this paper, we compute the correlation 
functions of the free fermion model by means of the algebraic Bethe ansatz method [1, 2, 3]. 
Our computation arc based on the recent progress on the Drinfeld twists. Working in the 
F-bases provided by the F- matrices (Drinfeld twists), the authors in [4, 5] managed to derive 
the determinant representations of the form factors and correlation functions of the XXX 
and XXZ models in the framework of algebraic Bethe ansatz. 

Recently we have constructed the Drinfeld twists for both the rational gl{m\n) and the 
quantum C/g(g'/(m|n)) supersymmetric models and resolved the hierarchy of their nested 
Bethe vectors in the F-basis [6, 7, 8]. These results serve as the basis of our computation in 
this paper of the correlation functions of the Uq{gl{l\l)) free fermion model. 

Correlation functions of the free fermion model based on the XXO spin chain (XY model 
[9]) with periodic boundary condition were studied in [10]- [15]. As is seen in section VI, by 
using the Jordan- Wigner transform, our Uq{gl{l\l)) free fermion model is equivalent to a 
twisted XXO model, and the one-point functions we obtained (see (V.5) and (V.7) below) 
give the m-point correlation functions of the twisted XXO model (see e.g. (VI. 6)). 

The present paper is organized as follows. In section II, we review the background 
of the Uq{gl{l\l)) model and its algebraic Bethe ansatz. In section III, we construct the 
Drinfeld twists for the model. In section IV, we obtain the determinant representation of the 
scalar products of the Uq{gl{l\l)) Bethe states. Then in section V, we compute correlation 
functions of the local fermionic operators of the model. We conclude the paper by offering 
some discussions in section VI. 

II Uq{gl{l\l)) free fermion model 
II. 1 The background of the model 

Let V be the 2-dimensional L^q(5f/(l|l))-module and R G End{V®V) the i?-matrix associated 
with this module. V is Z2-graded, and in the following we choose the FB grading for V, i.e. 
[1] = 1, [2] = 0. The i?-matrix depends on the difference of two spectral parameters Ui and 
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U2 associated with the two copies of V, and is, in the FB grading, given by 

/ ci2 \ 
D / \ u f \ ai2 6i2 

Rl2{Ul,U2) ^ Rl2{Ui- U2) = n h- H ' 

U Oi2 O-li U 

V 1/ 



(II.l) 



where 



ai2 = a{ux,U2) 
Cl2 = c(mi,U2) : 



sinh(wi — U2) 
sinh(wi — U2 + r]) 
sinh(Mi — U2 — f]) 



'12 



b'^{ui,U2) = - 



sinh(Mi — U2 + r]) 



(11.2) 



sinh('Ui — U2 + r]) 

with r] & C being the crossing parameter. One can easily check that the it!-matrix satisfies 
the unitary relation 

R21R12 = 1. (II.3) 

Here and throughout Rij = Rij{ui,Uj). The i?-matrix satisfies the graded Yang-Baxter 
equation (GYBE) 

R12R13R23 — R2zRi3Ri2- (II-4) 
In terms of the matrix elements defined by 



(11.5) 



the GYBE reads 



i',j',k' 

= E R{u2-u^^;i;R{u,-u^'R{u,-U2)^f{-Vf^'^^^^'^. (II.6) 

i',j',k' 

The quantum monodromy matrix T{u) of the free fermion model on a lattice of length 
N is defined as 



Tq{u) = Ron{u, Zn)RqN-i{u, Zn^i)_Roi{u, Zi), 



(11.7) 



where the index refers to the auxiliary space and {zi} are arbitrary inhomogeneous param- 
eters depending on site i. T{u) can be represented in the auxiliary space as the 2x2 matrix 
whose elements are operators acting on the quantum space V®^: 



Tn(u) 



A{u) B{u) 
C{u) D{u) 



(II.8) 



(0) 



By using the GYBE, one may prove that the monodromy matrix satisfies the GYBE 

Ruiu - v)T,(u)T2(v) = T2(v)T,(u)Ri2(u - v). (II.9) 

or in matrix form, 

= j:T{v)^Tiu)tRiu - („.io) 

i'd' 

Define the transfer matrix t{u) 

t{u) = stroTo{u), (11.11) 

where stro denotes the supertrace over the auxihary space. With the help of the GYBE, 
one may check that the transfer matrix satisfies the commutation relation [t{u),t{v)] — 0, 
ensuring the integrability of the system. The transfer matrix gives the Hamiltonian of the 

system: 



H = 



d\Rt{u) 
du 

N 



«=0 



1 

+ '"'-BS^^i))) . (11-12) 

where E'^^^ are generators, which act on the kth space, of the superalgebra Uq{gl{l\l)). 
Using the standard fermionic representation 

^11) = Cfe, Ell) = 4, Ell-) = 1 - Hfc, Ell) = rife, rife = c{ck, (11.13) 
the Hamiltonian can be rewritten as 

^ = ^i^ ^ (''^■''^^+^ + - 2coshr/(l - n,)) . (11.14) 

II. 2 Algebraic Bethe ansatz 

The transfer matrix (11.11) can be diagonalized by using the algebra Bethe ansatz. Define 
the Bethe state of the system 

11 

^n{v^,V2, . . . , ^;„) - n C(^0|0 >, (11.15) 
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where |0 > is the pseudo- vacuum, 



-»n(;) 



(11.16) 



(k) 



and the index (k) indicates the kth space. 

Applying the elements of the monodromy matrix (II. 8) to the pseudo- vacuum |0 > and 
its dual, we easily obtain 

B{u)\0>=0, <0|C(m) = 0, L'|0>=|0>, <0\D{u)=<0\, 

N N 

A{u)\0>^l[a{u,Zi)\0>, <0\A{u) ^l[a{u,Zi) <0\. (11.17) 

i=l i=l 

With the help of the GYBE (II. 9), we obtain the commutation relations between the 
elements of the monodromy matrix 



C{u)C{v) = 

D{u)D{v) = 

A{u)C{v) = 

D{u)C{v) = 

B{u)C{v) = 



-c{u,v)C{v)C{u), 

D{v)D{u), 

a[u,v) a{u,v) 

^ -C{v)D{u)-^-4^C{u)D{v), 



a{v, u) 



a(v, u) 



(11.18) 
(11.19) 
(11.20) 

(11.21) 



-C{v)B{u) + ^^^j^[D{v)A{u) - D{u)A{v)] 



-C{v)B{u) + 



a{u, v) 
b'^{u, v] 



[D{u)t{v) - D{v)t{u)]. 



(11.22) 



a{u, v) 

Thus applying the transfer matrix t{u) = D{u) — A{u) to the Bethe state and using the 
commutation relations repeatedly, we obtain the eigenvalues of t{u) as 



t(w)$iv = A(«, {^^fc})$ 



N 



a[Vk, u) f,^^ a[u, Vk) 



N 



providing Vk{k — 1,2, ... ,n) satisfying the Bethe ansatz equations (BAE) 

N 

n ^{'"k.zj) = 1. 

For late use, we define the state of the free fermion chain of length N 
\a\a2 ... Oat >= |ai >(i) |a2 >(2) • • • \aN >{n) 

and its dual 

\a1a2 ... Oat >'^=< aN\{N) < ctN-i\{N-i) • • • < =< QnCIn-i • • • 



(11.23) 



(11.24) 



(11.25) 



(11.26) 
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Ill Drinfeld twists of the model 



III.l Factorizing F-matrix and its inverse 

Following [4], we now introduce the notation jy, where a is any element of the permutation 
group Sn- We note that we may rewrite the GYBE as 

where To^23 = R03R02 and (T23 is the transposition of space labels (2,3). It follows that Ri j^ 
is a product of elementary i?- matrices [4, 6], corresponding to a decomposition of a into 
elementary transpositions. With the help of the GYBE, one may generalize (III.l) to a 
A^-fold tensor product of spaces 

Ri...n'^0,1...N = To^a{l...N)Ri,„Ni (111-2) 

where To^i,„n = Ron ■ ■ ■ Rqi- This implies the "decomposition" law 

Rl.^N — Ra'{l...N)Rl...Nj (III. 3) 

for a product of two elements in S^. Note that -R^/(i ^) satisfies the relation 

Ra'{l...N)'^0,a'(l...N) = To^a'a{l...N)Rl' (1...N)- (III-4) 

Let us write the elements of -Ri . jv as 

{rl.n);:^,''''''' , (ni.5) 

where the labels in the upper indices are permuted relative to the lower indices according to 
a. 

We proved in [6, 7, 8] that for the i?- matrix Ri j^, there exists a non-degenerate lower- 
diagonal F-matrix (the Drinfeld twist) satisfying the relation 

Fa{l...N){Za{l), • ■ ■ , Za(N))Ri...Ni^l^ • • • > ^n) — -^l...Ar(^l, ■ ■ ■ , Z^)- (III. 6) 

Explicitly, 

* N 

Fl,...^^= E E ]lPaif;'Sic,a,a,)Rl,,^, (111.7) 

aeSN <^a(l)—<^a(,N) j=l 
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where the sum Yl* is over all non-decreasing sequences of the labels aa{iy- 

aaii+i) > C(a{t), if cr{i+ 1) > a{i), 

Q!^(i+i) > acr{i), if + 1) < a{i) (III.8) 
and the c-number function S{c, a, a^) is given by 

Sic a, a,) = exp 1 1 (l - (-l)!"^^^)]) ln(l + c<,(feV(o) I • (III.9) 

( ^ l>k=l ) 

The inverse of the F-matrix is given by 



Fr\-K..Nll^^' (iii.io) 



i<j 



with 



and 



Aij = diag ((1 + Cij){l + Cji),aji, aij, 1) (III.ll) 



** N 
creSjv °ct(jv) j=l 



(III.12) 



where the sum is taken over all possible ccj which satisfies the following non-increasing 
constraints: 

aa{i+i) < Oia{i), if cr{i + 1) < 

< aa(i), if cr{i + l) > a{i). (III. 13) 

III. 2 Symmetric representation of the Bethe state 

The non-degeneracy of the F-matrix means that its column vectors form a complete basis, 
which is called the F-basis. By the procedure in [8], we find that in the F-basis, the simple 
generators of the superalgebra Uq{gl{l\l)) have the symmetric form: 



Pl2 _ 771 77ll2 I7I-1 

N 

= E 4) ^.-7^^ diag (26-" cosh 77,6-") (III.14) 

i=i 

£^21 771 77l21 771 — 1 

-c/ — Pi2...N-t^ i'i2...N 



Y: diag (e^{2aj, cosh r))-\ e%l) . (III.15) 



Similarly, the diagonal element D[u) of the monodromy matrix in the F-basis is given by 

D{u) = Fi2...iv£>(M)Fr2'..Ar = «)f=idiag(a„,-,l), (III.16) 

where aoj = a{u, zj). 

Then, the creation and annihilation operators C{u) and B{u) read, in the F-basis, 
C(u) = Fu...NC{u)Frl,j,^(q-'^^D{u)-D{u)^^)q-^l^'^^^^ 

N 

= E ^oi4) ^^^i ^i^g (2«o^- cosh V, 1)(,) , (111.17) 

i=l 

B{u) = Fu...nB{u)F^I j^ = gSti (E^^D - qDE^') 

N 

= - E Ki^i^ ®3^i diag (ao,(2a,i cosh ?7)-\ ajl) , (III.18) 
i=i ^•''^ 

where 6^- = -Zj), q = e^^ and /i = — — E"^^. 

Acting the F-matrix (III. 7) on the state (11.16), one sees that the pseudo- vacuum is 

invariant. Therefore in the F-basis, the Bethe state (11.15) becomes, 

n 

^N{Vl,V2,...,Vn) = Fi,„N'^N{vi,...,Vn) = l[C{Vn)\0 > . (III. 19) 

i=l 

Substituting (III.17) into (III.19), we obtain 

^n{vi, ...,Vn)^{2 cosh 77)"^ B-{vi, Vn\Zi„. . . , jFg) . . . Fg) |0 > , 

i\<...<in 

(III.20) 

where 

n n 

Bn{vi,...,Vn\zi^,...,Zi^) = ^ sign((7) [] ^^(^'fe, ^M*)) n a(^fe'^Mo) 

ceSn k=l l=k+l 

= deiB^{{vu],{z,]) (III.21) 
with B^{{vi\, {zj}) being a. n x n matrix with matrix elements 

B^^ = b^{v^, zp) n a{v„ zp). (III.22) 

7=1 

Similarly, acting B{un) . . . B{ui) on the dual pseudo-vacuum state, we have, 

- ( -1) n N 

<Q\B{un)...B{ur) = (-l)'^(2coshr7)^^ E 11 11 o.-\zk,z,^) 

i\<...<in 1=1 k=l,^ii 

xdetB^iM, {zij) < 0|^£) • • • Efl) ■ (111.23) 
8 



IV Determinant representation of the scalar product 
of the Bethe states 

In [2, 5], the authors gave the determinant representation of the scalar product of the Bethe 

state for the spin-1/2 XXZ model. In this section, wc derive the determinant representation 
of the scalar product of the Uq{gl{l\l)) Bethe states defined by 

Sn{{uj}, {vk}) =< 0|SK) . . . B{u,)C{v,) . . . CK) |0 > . (IV.l) 

The F-invariance of the pseudo- vacuum state |0 > and its dual state < 0| leads to 

Sn{{uj}, {vk}) -< 0\B{u„) . . . B{m)C{v,) . . . C{Vr,)\0 > . (IV.2) 

Following [5], we define 

, 111, . . . , Uniy • • • ) ^nj 

= <in,..., lm+l\B{Um) • • • B{Ul)C{v,) . . . C{Vn)\0 >, (IV.3) 

where ik {k = m + 1, . . . ,n), ordered as im+i < ■ ■ ■ < in, indicate the positions having state 
^ Q ^ , and other positions have state ^ 'i' ^ • One sees that when m = n, G*^"^ = Sn- 
Inserting a complete set and noticing (III. 18), (IV.3) becomes 

) ^1) • • • ) Uirm ■ ■ ■ j'^n) 

N 

— ^ ] ^ ^n; • • • ) l-^('^m) |^m+l) ■ ■ ■ > im-\-pi ji ^m+p+l; • • • ) ^ 

). (IV.4) 

In view of (III. 18), we have 

< ijj, . . . , im+l \B(^Um^ linfi+l, ■ ■ ■ ■> im+pi ji ^m+p+li ■ ■ ■ tin ^ 

N n 

= -{2coshriP^""^^-i-irb-'iu^,Zj)l[a-\zk,Zj) H a{um,z,,). (IV.5) 

k^j l=m+l 

With the help of (III.20), we obtain 

n 

G^^\{vk},ii,...,in) =< ^n,•••,^l| n<^K)|0 > 

k=l 

= (2coshr7)'^detB-(K},{;^J). (IV.6) 
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We now compute G'^^^ by using the recursion relation (IV. 4). Substituting (IV. 5) and 
(IV. 6) into (IV.4), we obtain 



G^^\{vk},ui,i2, ■■■,in) 

N 

= H < ^n, • • • ,^2|5(Ml)|i2, • • • , ip+l'i' W2, • • • > 

3^12, ■■■,in 

xG^^\{vk},i2, ■ ■ ■ ,ip+i,j,ip+2, • • • ,4) 

(n-l)(n-2) ^ N n 

= -(2cosh77) 2 J2 {-iyb+{ui,Zj)'[[a~^{zk,Zj)Yla{ui,Zi^) 

j=ii2,...,in kjtj 1 = 2 

xdetB~{{vk}, Zi^, Zi^^„Zj, Zi^^^, ZiJ, {k = l,---,n). (IV.7) 

Let Vk{k — 1, . . . , n) label the row and zi{l — i2, . . . , j, . . . , in) label the column of the 
matrix B". Prom (IV. 6), one sees that the column indices in (IV.7) satisfy the sequence 
i2 < . . . < j < ■ ■ ■ < in- Therefore, moving the column j in the matrix B~ to the first 
column, we have 

G^^\W},Ui,i2, ...,in) 

(n-l)(n-2) ^ N n 

= -(2cosh?7) 2 J2 b-^iui, Zj)'[la~^{zk,Zj)Y[a{ul,Zi^) 

jf^h,...,i„ k^j 1=2 

xdeiB{{vk},Zj,Zi^,...,ZiJ 
= -(2 coshry)^''^^ det(i3-)(^)({^fc}, u^, z^,, . . . , z^J, (IV.8) 

where the matrix {B~Y^\{vk} ■, ui, Zi^, . . . , Zi^) is given by 

(B-^)« = a(u,,z,^)B-f, for/3>2, (IV.9) 

N 0-1 N 

i^aif^ - E b^{uuz,)b-{v^,zj)l[a{v^,z,) [] <^-\^k,Zj). (IV.IO) 

j^i2,--;in 7=1 k=l,^j 

Using the properties of determinant, one finds that if j = i2, • • • , in, the corresponding terms 
in (IV.IO) contribute zero to the determinant. Thus, we may rewrite (IV.IO) as 

^ V e"^--°sinh^r? sinh(^;^ - ^,) 



sinh(Mi - Zj + T]) sinh(i;Q, - Zj + rj) sinh(u^ - Zj + rj) 
^ sinh(zfc — Zj + rj) 



X 



n 



smh{zk - Zj) 

= e«V(«i)- (IV.ll) 
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Thanks to the Bethe ansatz equation (11.24), we may construct the function 



e 



sinh(vc, - Up) sinh(w^ - up) \ smh{up - Zk + rj) j ' 

Comparing f{ui) in (IV. 11) with g{ui) in (IV. 12), one finds that as functions of Mi, they 
have the same residues at the simple pole ui = zj — r] mod(i7r), and that when ui — >■ cxo, they 
are bounded. Moreover, one may prove that the residues of g{ui) at ui = {ly = 1, . . . , a) 
are zero because Vi, are solutions of the Bethe ansatz equation (11.24). Therefore, we have 

iB-^)^'^ = -^-}^^I[a-\u„v,) (l-fi a{u^,z,)] . (IV.13) 

Then, by using the function G^^^ , G^^^ and the intermediate function (IV.4) repeatedly, 
we obtain G^'^^ as 

, 111 J • • • . '^^'rm- '^'rn-\-l: ■ • • 5 ^nj 

ni n — 1 ) — ni( 2 (! — fn — 1) _ — _ 

= (-l)™(2cosh77) ^ n ^ i^k,Uj) 

l<j <k<Tn 

Xdet{B-)^"'\{Vk},Ui, . . .,Um,tm+U • • • , ^n) (IV.14) 

with the matrix elements 



fc=l 

i^ap)^""^ = M-p, for/3<m. (IV.15) 



(IV.14) can be proved by induction. Firstly from (IV.8), (IV.9) and (IV.13), (IV.14) is true 
for m — 1. Assume (IV.14) for G^"^~^\ Let us show (IV.14) for general m. Substituting 
Q{m-i) g^^^ (IV.5) into intermediate function (IV.4), we have 

G'(-)(K} , Itl, . . . , UfYi, Irn+lj • • • ) ^nj 
N 

— ^ ] ^ ini ■ ■ ■ 1 im+1 \ B{Urri) \ im+l ) ■ ■ ■ i ^m+pi ji ^m+p+l; • • • ) ^ 

XG^ ^ (^{Vk^j Ui, . . . , Ufn—i, ijn-\-\, . . . , im+pj ji • • • ) in) 

N N n 

= -(2cosh?7)-("-™) ^ b-^{um,zj)l[a-\zk,zj) [] a{um, ZiJ 

j^im+i,—,in k^j l=m+l 

XG'("^-1)({V^}, Ml, ... , Um-l,j, im+1, ■■■,in) 
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. rt(n—l)—m(2n—m—l) — r i 

(-ir(2coshr7)^ ^ n a-\uk,Uj) 

l<j<k<m~l 



where the matrix elements B'^^^ are given by 



xdetB'^'^'{{vk},ui . . .,Um,im+l, ■ ■ ■,in), (IV.16) 

m) 
m 

^'i? = n ^i'^'^' for 13 >m, 

k=l 

By = M-^ for/3<m, 

m— 1 o— 1 

^Im = n«("i'^J-) h-^{u^,Zj)h-{va,Zj)^a{v^,Zj) 

Af 

X n (^~\^k.z^). (IV.17) 
Thus, by the procedure leading to we can prove 

m—l 

^'2= n «.)A^-^. (IV.18) 

Then one sees that B'^^j^ = Therefore we have proved that (IV. 14) holds for all m. 

When m = n, we obtain the scalar product Sn{{ui}, {vj}), 

n 

Sn{{u^},{vJ}) = (-1)" n«"'K'«OdetM-({^;,},M), (IV.19) 

k>l 

where the matrix elements of are given by 

M% ^^-^I^^Y[a-\up,v,) fl - n a{u„z,)] . (IV.20) 

By using the expression of the eigenvalues of the system (11.23), the scalar product (IV.19) 
can be rewritten as 

n 

'5n(M,{^.}) = (-l)"n«^'K>«Odet^-({t;,},^.}) (IV.21) 

k>l 

with the matrix being 

= e±(''«-"'^) sinh(«^ - v^) n a(^;^, u^) H a-\u^, ^^^dj±^{v^_ ^^^22) 

Remcirk: In the derivation of (IV.19), the parameters Vi in the state C{y^ . . .C'(v„)|0 > 
are required to satisfy the BAE (11.24). However, the parameters {j — 1, . . . ,n) in the 
dual state < Q\B{un) . . .B{ui) do not need to satisfy the BAE. 
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On the other hand, if we compute the scalar product by starting from the dual state 
< 0\B{vn) . . . B{vi), then by using the same procedure, we have 

Sn{{vi},{uj}) = <0|BK)...B(^;i)C'M...C'K)|0> 

n 

= {-irl[a-\uk,ui)detM^i{v,},{uj}). (IV.23) 

k>l 

In the above equation, we have assumed that {vi} satisfy the BAE. 



Noticing the BAE (11.24), one sees that the scalar product Sn{{ui} , {vj}) = if both 
parameter sets {ui} and {vj} {{vj} ^ {ttj} z, j = 1, . . . , n) in (IV. 19) and (IV.23) satisfy the 
BAE. 

Let Ua ^ Va {a = 1, ... ,n) in (IV. 19), we obtain the Gaudin formula for the norm of 
the Ug{gl{l\l)) Bethe state. 



Sn = Sn{{Vj},{Vk})=<0\B{v,,)...Bivi)Civi)...CiVn)\0> 



k>j 

= (-l)"sinh"r/n 



= (-l)"sinh'"r7n 



sinh^(t;fc — 


Vj + T]) 


sinh^(ffc 


-Vj) 


sinh^(t'fc — 


Vj + r]) 


sinh^(fjk 




sinh^(i;fc - 


- Vj + 1]) 



n 



.0=1 



'-j^ _ -Q sinh(M„ - zi) 



sinh(M„ -zi + 7])J\ 



TT In '^^'^^(^" ~ ^i) 

iJi dua Vi=i sinh(Ma - Zi + T})^ 



n N 

HE 



j^^- sinh {vk - Vj) sinh(^;a - zi) sinh(^;c. - zi + r]) 

where we have used the BAE (11.24). 



, (IV.24) 



V Correlation functions 

Having obtained the scalar product and the norm, we are now in the position to compute 
the k-point correlation functions of the model. In general, a k-point correlation function is 
defined by 



'-'^ =< 0\B{u^) . . . B{m)el . . . e^Civ,) . . . C(^„)|0 >, 



(V.l) 



where e^^ stand for the local fermion operators q^. , cj^ or rij^. , and the lower indices ij indicate 
the positions of the fermion operators. 
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The authors in [16] proved that the local spin and field operators of the fundamental 

graded models can be represented in terms of monodromy matrix. Specializing to the current 
system, we obtain 

j-i N 
4 = Ili-M^k) + D{zk)) ■ B{zj) ■ H {-A{zk) + D{zk)), (V.2) 

k=l k=j+l 
i-1 N 

Ci = ]li-M^k) + D{z,))-C{zj)- II {-A{z,)+D{zk)), (V.3) 

k=l k=j+l 
7-1 N 

^3 = Y[{-A{zk) + D{zk))-D{z^)- {-A{zk) + D{zk)). (V.4) 

fe=l k=j+l 

V.l One point functions 

In this subsection, we compute the one point functions for the local operators cj^^, Cm and 
nm, respectively. 

We first calculate cj^^ . Noticing that the Bethe state and its dual are eigenstates of the 
transfer matrix under the constraint of the BAE, we have, from (V.2), 

({%}>^m,K}) 

= < 0|BK) . . . B{ur)clC{v^) . . . CK+i)|0 > 

= ^m-l{{Uj})(t>-l{{Vk}) < 0|BK) . . . B{ui)B{Zrn)C{v{) . . . C{Vn+l)\0 > 
= cl>rn-l{{uj})<j);;,\{Vk})<0\B{Un)...B{u,)B{Zrn)C{vi)...C{^^^^ 

= ...,1*1, Zm, {Vj}) 

n n 

= (-l)"+Vm-i({Mi})C^(K}) n a'\uk,Uj)Y[a'\ui,Zm) 

k>j 1=1 

XdetM~{{Vj}, Zm,Ui,..., Vn), (V.5) 

where 0j({wj}) = 111=1 nr=i <^{ui,Uk)- As mentioned in the remark of the previous section, 
F~ = if the parameter set {-Uj} (i = 1, . . . ,n) is not a subset of {vj} {j = 1, . . . , n + 1). 
When {ui} C {vj}, (V.5) can be simplified to a simple function. For example, if tij — '^i+i 
{i — 1, . . . ,n), the one point function F~ becomes 

Vi, . . . ,fn+l) 

, ^ ^+i 0„.-i({»,})r-('--^---^smli^"+^7y smli^(rA:-r,+7/) '^ smli(r,-.,„ + r^) 

(t>m{{vk}) sinh(i;i - Zm) k>j^2 sinh^(^;fc - Vj) jj^ sinh(t;j- - Zm) 

^/j±. ^inh(.,-.,+,) -^^^ 1 

j-J-^ smh{vj - Vi) smh{va - Zi) smh(i;a - Zi + r]) 
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Similarly, when {ui} C {vj}, we obtain the one-point function involving the operator Cm'- 

= < 0\B{vn+i) . . . B{v,)c^C{m) . . . C{un)\0 > 

= (j)m-l{{Vj})(l);;,^{{Uk})Sn+l{{Vj}, Zm, Ui, . . . , Un) 

n n 

fc>i i=\ 

xdetA^+({vj}, -ui, . . . , Vn). (V.7) 

is non-vanishing if {wj} C \yj\- When {wj} C {wj}, (V.7) can also be simphfied to a 
simple function. In the case ttj = Wj+i (z = 1, . . . , n), the one point function becomes 

F+(v„+i, . . . , t;2, ^„j, vi, . . . , 

. e(^i-^-) sinh^"+^ r] sinh^(^fc - + rj) sinh(^,- - -g^ + 77) 

<^m{{uky) sinh(i;i - z„) ^,^^:L2 sinh^(vfe - Vj) fj^ sinh(i;j- - Zm) 

^/+^ smh(.,-.,+,) -|+^^ ^ 

^^^^2 sinh(vj - vi) sinh(t;o - Zi) smh{va - Zi + r]) 

The one-point function involving the operator rim is defined by 

F:-{{u,}, Zm, {Vk}) =< 0|SK) . . . B{m)nmC{vi) . . . C{vr^+,)\0 > . (V.9) 
Substituting (V.4) into the above equation and considering the BAE, we have 

Zm, {vk}) =< 0\B(un) . . . B(m)nmC(vr) . . . CK)|0 > 

(l)m-l{{Uj}) 



0m-l(W}) 

With the help of (11.21), we see 



< 0\B{Un) . . . B{ui)D{zm)C{vi) . . . C{Vn)\0 > . (V.IO) 



D{Zm)C{v,)...C{Vr,)\0> 

n 

l[a-\vk,Zm)C{v^)...C{vr,)\0> 



k=l 

J-1 



XC{VI) . . . C{vj^i)CiZm)C{vj+i) . . . C{Vn)\0 > . 



(V.ll) 
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Therefore, substituting (V.ll) into (V.IO), we obtain 

0m-l(K}) 

-2^—( tIIt V 11 " {Vk,Vj)Sn{{Ui},Vi,...,Vj_i,Zm,Vj+i,...,Vn) 

j=l (^K^ji^m) C[Vi, Zjn) /.^i -^j 

xdet [m+({m,}, {vj}) - ^f{{ui}, {vj}, zm)] , (V.12) 

where jV is a rank-one matrix with the following matrix elements 

xr ff \ f \ \ e""-^'3sinh^r7 sinh(z^ - + r/) 

smh{ua - Zm) smh{vf3 -Zm + V) t=i smh(^„ - Ui) 

In the above derivation, we have used the following property of determinant: If A is an 

arbitrary n x n matrix and B is a rank-one n x n matrix, then the determinant of .4. + i3 is 
given by 

n 

det (A + B)^ detA + det^(^) , (V. 14) 

i=l 

where 

A% = Aap for (3^1, 
A^^- = B ■ 



V.2 Correlation function of two adjacent operators 

In the subsection, we compute the correlation function of two adjacent operators Cm and 

Cm+i defined by 

1 Zm-i Zm+\i 

{vj}) =< 0|SK) . . . B{u^)cmcl+^C{vr) . . . C(t;„) |0 > . (V.15) 

Substituting (V.3) and (V.2) into the above definition and considering the fact nfeLi^(^fc) = 
1, we have 

= <0\B{Un)...B{u^)C{Zrn)B{Zrn+l)C{v^)...C{Vn)\0> . (V.16) 
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By using the commutation relation (11.22), we obtain 

B{zr,+,)C{v,) . . . CK) |0 >= i-irCiv,) . . . C{vn)B{zm+i)\0 > 

j^l a[Zm+l,Vj) 

+ J2{-iy^-^j^^^^^^C{v,) . . . C{v,_,)t{^^+i)D{v,)C{v,+,)C{vr.)\0 >, (V.17) 

j^l a[Zm+l,Vj) 

where i(u) = Fi,,,Nt{u)F^^j^. On the rhs of the above equation, one easily finds that the 
first term is zero. Using the BAE, one may check that the second term also equals to zero. 
Therefore, only the third term survives on the rhs of the above equation and we have 

B{Zm+l)C{Vi)...C{Vn)\0> 



n 



xC{v,) . . .C{vj_,)C{v,+,)C{v^)\0 > 

j=l "V^m+l: Uj) k=j+l 

" b-{vi,Vj) '-pi c{v„„vi) -A- ^1 

l=j+l '^y^l' '^0) m=j+l ^y^rn-, Uj) ^^j^^^^i 

xC{vi) . . . C{vj.i)C{vj+,) . . . C{vi^i)C{v,)C{vi+,) . . . C{vn)\0 > 

n 

= Y.M,C{v,) . . .C(^;,_i)C(^;,+i)CK)|0 > 

a n 

+ E E MJ,lC{v^) . . . C{vj.,)C{vj+,) . . . C{vi^,)C{v,)C{vi+,) . . . C^) |0 > . 
j=i i=j+i 

(V.18) 

Substituting (V.18) into (V.16), we obtain two-point correlation function 

4>m-l{{Ui}) 



(pm+l{{Vj}) 



J2 MjSn{{Ui}, Zm,Vi,..., Vj-i, Vj+i, Vn) 
3=1 



+ E E Mj^iSn{{Ui}, Zm,Vi,..., Vj_i,Vj+i, . . . , Vi_i, Vj, Vi+i, . . . , W„) 
j=l l=j+l 



. (V.19) 
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VI Discussion 



In this paper, with the help of the factorizing F-matrix (F-basis), we have obtained the de- 
terminant representations of the scalar products and correlation functions of the Uq{gl{l\l)) 
free fermion model. 

In [10]- [15], the authors studied the correlation functions of the free fermion model based 
on the finite XXO spin chain (XY model [9]) with periodic boundary condition 

N 

Hxxo = E (^I^I+i + ^l^l+i + ha]) , (VI. 1) 

j=i 

where {e — x,y,z) are the Pauli matrices and h is an external classical magnetic field. 
The equivalence between the free fermion model and the XXO model can be proved by using 
the Jordan- Wigner transform 

Cfe = exp[i7r(5jt_i](7jf , (VI.2) 
4 = exp[iTrQk-i], (VI.3) 

where — |((T^ ± cr^), Qk — J2'j=i i(l ~ <^fe)- Because of the periodic boundary condition 
of the finite XXO spin chain, we have 

4+1 = 4- (VI.4) 
Substituting the Jordan- Wigner transforms into the above relation, we obtain 

cn+1 = exp[z7r(5Ar]ci, cj^+i = cl expfiTrQjvj- (VI.5) 

Thus, comparing the above boundary condition with that of the Uq{gl{l\l)) free fermion 
model (11.14), we find that the free fermion model arising from the XXO model has a twisted 
boundary condition which depends on the operator = Y,i^=i ^l- 

On the other hand, by means of the Jordan- Wigner transform, the Uq{gl{l\l)) free fermion 
model is equivalent to a twisted XXO model, and the one-point correlation functions (V.5) 
and (V.7) give rise to the m-point correlation functions of the twisted XXO model. For 
example: substituting (VI.3) into (V.5), we obtain 

-^n ({«j}>^m,{^^fc}) 

= < 0[S(m„) . . . B{u^)clC{v^) . . . C(^„+i)[0 > 

= < Q\B{un) . . . B{ui)al . . . <_ia-C(^;i) . . . C{vn+i)\Q > • (VI.6) 
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